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Abstract: 

We analyze the computational power and limitations of the recently proposed 'quantum adiabatic evolution algorithm' . 

1 Introduction 

Quantum computation is a revolutionary idea that has fundamentally transformed our notion of feasible computation. 
The most dramatic example of the power of quantum algorithms was exhibited in Shor's celebrated quantum algo- 
rithms for factoring and discrete log [p"3|], Grover's quantum search algorithm [|l0| gives a quadratic speedup for a 
much wider class of computational problems. Despite numerous attempts in the last few years, it has proved to be a 
difficult challenge to design new quantum algorithms. Recently, Farhi et al. || proposed a novel paradigm for the 
design of quantum algorithms — via quantum adiabatic evolution. This paradigm bears some resemblance to simu- 
lated annealing, in the sense that the algorithm starts from an initial disordered state, and homes in on a solution (by 
what could be described as quantum local search) as a parameter V is smoothly varied from to 1. The challenge 
lies in showing that the process still converges to the desired solution with non-negligible probability if this transition 
is made in polynomial time. In [Q, 0], this paradigm was applied to the Exact Cover problem (which has a close 
connection to the 3SAT problem), and using computer simulations it was shown that the algorithm works efficiently 
on small randomly chosen instances of this problem. 

In the first part of the article, we discuss the quantum adiabatic theorem and explain the quantum adiabatic approach 
to computation. Next, we clarify the connection between the continuous time evolution of adiabatic computing and 
the quantum circuit model with its discretized time. We do this by describing a way of efficiently simulating quantum 
adiabatic algorithms with a network of standard quantum gates. After this exposition, we explore three questions about 
quantum adiabatic evolution algorithms. 

Can we apply the exponential lower bounds for quantum search [Q] to conclude that the adiabatic quantum algo- 
rithm for 3SAT must take exponential time? More concretely, at a high level of abstraction, the adiabatic quantum 
algorithm for 3SAT may be viewed as some quantum process that gets information about the 3SAT instance only 
by (quantum) queries of the following type: given a truth assignment, how many clauses of the formula $ are not 
satisfied? We prove that there is a (classical) polynomial time algorithm that can reconstruct the 3CNF formula $ by 
making polynomially many queries of this type. It is somewhat surprising that this question does not appear to have 
been studied in the context of relativization results for NP. In our context, it rules out any query complexity based 
(quantum) lower bound for the adiabatic quantum solution of 3SAT. 

Is adiabatic quantum computing really quantum? We give an example of an adiabatic quantum algorithm for 
searching that matches the optimal quadratic speedup obtained by Grover's search algorithm. This example demon- 
strates that the 'quantum local search', which is implicit in the adiabatic evolution, is truly non-classical in nature from 
a computational viewpoint. 

Finally, we give a simple example of a computational problem on which the adiabatic quantum algorithm provably 
takes exponential time. Although the problem is easy to solve classically, it is designed to be difficult for algorithms 
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based on local search: its global optimum lies in a narrow basin, while there is a local optimum with a much larger 
basin. Let / be a function on the n-bit strings, where f(x) depends only on w(x), the Hamming weight of x. The 
problem is to find an x that minimizes f(x). (Obviously, it is straightforward to solve this class of problems in n + 1 
steps.) Consider functions / such that for w(x) < (h + s)n, f(x) — w(x), and which decreases for w(x) > (| + e)n 
to the global minimum /(!" ) = — 1. We prove that for such instances, the adiabatic quantum algorithm requires an 
exponential slowdown in n. We do this by showing that the gap between the minimum and second eigenvalue of the 
Hamiltonian of the system is exponentially small. In an upcoming paper [ra], we generalize these techniques to show 
a similar exponential slowdown for 3 SAT. 

2 The Quantum Adiabatic Theorem 

The Hamiltonian of a physical system gives a complete specification of the time evolution of this system. At a given 
time t, let ip(t) denote the state of the system under the influence of the Hamiltonian H(t). The differential equation 
that describes the time evolution is the well-known Schrodinger equation: 

m£m)) = H(t)m)), 

where h is Planck's constant h w 6.63 x 10~ 34 Joule-second, divided by 2ir. A Hamiltonian is described by a Her- 
mitian matrix, whose eigenvectors represent the eigenstates of the system. The corresponding eigenvalues refer to the 
different energies of the eigenstates. The state (eigenvector) with the lowest energy (eigenvalue) is called the 'ground 
state' of the system. The Schrodinger equation can also be described with reference to the unitary transformation U 
that is defined by the Hamiltonian H (t) (from now on we work with h = 1): 

f t U{t) = -iH(t)U(t), 

with the initial condition U (0) = /. We say that the Hamiltonian evolution from H (0) to H(T) induces the unitary 
transformation U(T). The evolution of a system with a time-independent Hamiltonian H is easily expressed by the 
exponential U (T) = c~ lTH . Finding the (approximate) solutions for Hamiltonians that vary in time is one of the core 
tasks in quantum physics. One of the most important cases of such a time-dependent case is described by the adiabatic 
evolution of an isolated quantum mechanical system. 

The quantum adiabatic theorem states that a physical system that is initially in its ground state, tends to stay in this 
lowest energy state, provided that the Hamiltonian of the system is changed 'slowly enough'. 

The quantitative version of the adiabatic theorem gives the following specific upper bound on the slowdown that is 
required for the adiabatic evolution of the ground state. (See for example JT^ ] for more details on this.) Parameterize 
the time-dependent Hamiltonian by H(s) for < s < 1 and its ground state by <f>(s). Our goal is thus to gradually 
transform the applied Hamiltonian from H(0) to H(l) such that the initial state ^>(0) = (f)(0) evolves to a close 
approximation ip(l) ~ 0(1) of the ground state of H(l). We introduce a delay factor r(s), which determines the rate 
at which the Hamiltonian is modified as a function of s. Now the Schrodinger equation in s equals 

= -iT(8)H(8)M8)). 

The crucial quantity for this transformation is the gap between the two smallest eigenvalues of H(s), which we denote 
by g(s). It can be shown that a delay schedule r with 

T S > - 

is 'sufficiently slow' for the adiabatic evolution from <fi(0) to 4>(1). As a result, the total delay of this process will be 
of the order J s=(j r(s)ds. For most Hamiltonians it is too difficult to determine the gap g(s) for every s. If this is the 
case, we can also look at the minimum gap g m [ n := min s g(s) and the maximum A max := max s ||| j-£f(s)|||2, and 
obtain the adiabatic evolution with the constant delay factor r(s) = t c 6 0( f ax ). 



3 Adiabatic Quantum Computation 



Adiabatic quantum computation, as proposed by Farhi et al.[g], works as follows. At time t = 0, the quantum 
mechanical system is described by a Hamiltonian Hq, whose eigenstates are easy to compute. Next, this system is 
slowly transformed to its final Hamiltonian Hf, for which the ground state is the solution to a specific minimization 
problem /. We do this is by letting the energies \ z of the eigenstates z of Hf correspond with the function that we try 
to minimize. Hence, if this function / has domain {0, 1}™, then the final Hamiltonian is defined by 

H f ■= E /(*H*X*I- 

Z£{0,1}™ 

We will assume throughout this paper that / : {0, 1}* — > K is computable in polynomial time, and that f{x) is 
bounded by a polynomial in \x\. 

The choice of the initial Hamiltonian Hq is independent of the solution of the problem, and will be such that Hq is 
not diagonal in the computational z -basis. Specifically, we consider the 'Hadamard basis' with the bit values 

|6>:=^(|0> + |1» and |1> := J=(|0> - 

For a binary string z E {0, 1}™, let \z) denote the state which would be written as \z) in this basis. (The uni- 
tary mapping between these two representations is provided by the n-fold Hadamard matrix: W® n |2:) = \z) and 
W® n \z) = \z).) 

A simple starting Hamiltonian that fulfills the above requirements is 

H := £ h(z)-\z)(z\, 
ze{o,i}™ 

with h(Q n ) = and h(z) > 1 for all other z ^ n , such that the ground state with zero energy of Hq is the uniform 
superposition |0 • • • 0) = Yl z l z )- Having defined the initial and final conditions of our system, we will now 
describe the time-evolution. 

Following the proposal by Farhi et al. in [§, 0], we can define the time dependent Hamiltonian H(t) as the linear 
combination of the starting and the final Hamiltonian: 

H(t) := (l-±)Ho + ±H fl 

with < t < T, and T the crucial delay factor of the Hq — > Hf transition. 

By the adiabatic theorem we know that this system will map the initial ground state | -0 (0)) = |0") to the global 
minimum of the function /, provided that we pick T large enough. In the previous section we mentioned that T G 
0(A max g m ^ n ) is a sufficient upper bound on this delay. Without any further knowledge about the specific Hamiltonian 
H{t) — which involves detailed knowledge about the function /, this is also a lower bound for a reliable adiabatic 
evolution from Hq to Hf. Because ||| j^H(s) f% is polynomial in n (as long as / S poly(n)), we will ignore this factor 
and focus mostly on the T » g~ in requirement for the delay of the adiabatic quantum computation. 

4 Approximating the Adiabatic Evolution 

In this section we explain how the continuous time evolution from Hq to Hf can be approximated by a quantum circuit 
of size poly (nT). Our goal is to demonstrate the ingredients of the polynomial upper bound, and we do not try to 
optimize to get the most efficient simulation. 

The approximation is established in two steps. First, we discretize the evolution from Hq to Hf by a finite 
sequence of Hamiltonians H[, H' 2l . . . that gives rise to the same overall behavior. Second, we show how at any 
moment the combined Hamiltonian H'j = (1 — s)Hq + sHf can be approximated by interleaving two simple unitary 
transformations. 

To express the error of our approximation, we use the £2 induced operator norm "||| o || 2 ": 

|M| 2 := max |Mz| 2 . 

II 3= II 2 = 1 

The next lemma compares two Hamiltonians H(t) and H'(t) and their respective unitary transformations U (T) and 

U'(T). 



Lemma 1 Let H(t) and H'(t) be two time-dependent Hamiltonians for < t < T, and let U(T) and U'(T) be 
the respective unitary evolutions that they induce. If the difference between the Hamiltonians is limited by \H(t) — 
H 1 {t)\2 < 5 for every t, then the distance between the induced transformations is bounded by \U{T) — E/'(T)||| 2 < 
V2T5- 

Proof: Let ip(t) and ip'(t) be the two state trajectories of the two Hamiltonians H and H' with initially ?/>(0) = ip'(0). 
Then, for the inner product between the two states (with initially (tp' (0)\i/j(0)) = 1), we have 



www)) = -iWimH^-H'mm)- 



Because at any moment t we have || |-^(*)> ||a = WW^h = 1 and || jET (t) — ^Z"'(t) ||| 2 < 6, we see that at t = T the lower 
bound \(ip' (T)\ip(T))\ > 1 - TS holds. This confirms that for every vector ijj we have \\U(T)\ip) - U'(T)\tp)\\ 2 < 
V2TS- ■ 

This lemma tells us how we can deviate from the ideal Hamiltonian Hit) := (1 — i;)Ho + ^Hf, without 
introducing too big of an error to the induced evolution. As mentioned above, we will approximate the continuous 
H(0) — ► H(T) trajectory by a sequence of r Hamiltonians H[, . . . , H' r , each of which applied for a duration of — . 
This yields the unitary evolution U'(T), defined by 

U (1 ) ■= e r ' r . . . e ». r i i ( 

with for any 1 < j < r the Hamiltonian H'j := H{^-) = (1 — ^)Ho + (^)Hf. If we view H' as a time-dependent 
Hamiltonian H'(t) := H j[t) withj(t) := ], then we have the bound \\H{t)- H' (i)|| 2 < i||if/-ifo|||2 6 0{n d /r) 
for all t. By the previous lemma we thus have the bound \U{T) - f7'(T)|| 2 € 0{yjTn d /r). 

The second part of our approximation deals with the problem of implementing the unitary transformations U'j 
defined by 

jji ■- e -i£(i-i)H -i£(i)H f _ 

with elementary operations. 

The Campbell-Baker-Hausdorff theorem[|]| tells us how well we can approximate 'parallel Hamiltonians' by con- 
secutive ones: |||e' 4+B — e e B ||2 € 0(|| ABIJ2). Hence in our case, by defining 

U" := e -i$V--i)Bo . e -i?(f)H/ ) 

we get the approximation \\U^ - C/j'|| 2 € 0(^-|||Fo-ff/|||2)- This leads to \\U'(T) - f7"(T)||| 2 € 0(n d+1 T 2 /r), and 
hence also for the original transformation: \\U(T) - U" (T)\\ 2 E 0{n d+1 T 2 /r). 

Because Ho = J2 Z i s diagonal in the Hadamard basis {0, 1}™, and Hf = f( z )\ z )( z \ i s diagonal in 

the computational bases, we can implement the above U" as 

\j>> = w® n . f Q j . w® n . Ffj t 
with W® n the n-fold Hadamard transform, and F and Ff the appropriate phase changing operations: 

Foj\z) ■= e-^-m*) lz): 
F fJ \z) := e-^m*)\ z ). 
Because h(z) and f(z) are easy to compute, so are Fo and Ff. We have thus obtained the following theorem. 

Theorem 1 Let Hq and Hf be the initial and final Hamiltonians used in an adiabatic computation, with the function 
f G 0(n d ). Then, the unitary transformation U(T) induced by the time-dependent Hamiltonian H(t) := (1— A)flo + 
ipHf can be approximated by r consecutive unitary transformations U[' , . . . , U" with r £ 0(T 2 n d+1 ). Furthermore, 
each I!" has the form W® n FoW® n Ff and can thus be efficiently implemented in poly(nT) time. 

It is interesting to note that the W® n FoW® n Ff transformation has the same form as the 'Grover iteration' of the 
standard quantum search algorithmjl0[]. More recently, we also learned that the work of Hogg on quantum search 
heuristics [O] describes essentially the same algorithm as the adiabatic approach to minimization. 



5 Quantum Adiabatic Searching 



One question that should be asked first is if adiabatic quantum computing is truly quantum computing. In this section 
we answer this question affirmatively by reproducing the quadratic speed-up of Lov Graver's search algorithm. 

For the search problem, the function / : {0, 1}™ — > M takes on value 1 on all strings except the solution u E {0, 1}™ 
for which /(it) = 0. Thus the final Hamiltonian for the adiabatic algorithm, H u , will have eigenstates \z) with 
eigenvalue 1, with the exception of the unknown solution u 6 {0, 1}™, which has eigenvalue 0: 

H u ■= £ \z){z\. 

*6{0,1}"\M 

The initial Hamiltonian is defined similarly, except that it is diagonal in the Hadamard basis, and has ground state 

|6"): 

H := ]T \z)(z\- 

*e{o,i}«\{0"} 

With these initial and final conditions one can easily show that for the resulting time-dependent Hamiltonian 

Hit) := (1 — j<)Hq + y-ffn, 

the gap between the two smallest eigenvalues as a function of s := £ is expressed by 



9(>) _ ^(^-IX,*-^ m 
This gap reaches its minimum at t = ^ when it equals -y=- At first sight, this would lead to the conclusion that the 



necessary delay factor T = ^(5 m f n ) is linear in N = 2". However, by using our knowledge of the gap function g(^) 
we can significantly reduce the running time to 0{sfN). 

For example, regardless of the solution u, we know that the transition from H(0) to will have a minimal 

gap that is significantly bigger than -j=. The necessary delay factor that we use for this first part of our transformation 

Hq — > H u , can therefore be much smaller than N, In general at any moment s = ^, Equation [j] tells us the size of 
the gap g(s), and hence the delay factor that suffices at that moment. This means that we can employ a varying delay 
factor g(s)~ 2 , without destroying the desired adiabatic properties of the evolution Ho — > H u . In sum, this approach 
leads to a total delay factor of 

ds f 1 2" 

-ds 



o 9{s) 2 Js=o 2" + 4( 2"-l) ( S 2 - s) 

2" • arctan(V2 n - 1) 
V2^T ' 



As a function of N = 2", this gives a time complexity 0(y2 n ) = O(vN), which coincides with the well-known 
square root speed-up of quantum searching. (See the article by Farhi and GutmannlQ for another example of a 
'continuous time algorithm' for quantum searching.) 



6 Query Bounds for the 3SAT Problem 

The adiabatic quantum algorithms of [0, [| work on 3SAT as follows: on input a formula $ = C\ A • • • A Cm (where 
the Ci are clauses in variables xi, . . . , x n ), the only way the quantum algorithm gathers information about $ is by 
queries which ask, for a given truth assignment b (in general a superposition of assignments), how many of the M 
clauses b does not satisfy. A natural approach to establishing a lower bound on the running time of the adiabatic 
quantum algorithm is to show that any quantum algorithm must make a large number of such queries to solve the 
problem. This is the approach that leads to the exponential lower bound for unstructured search [Q] (there the query 
asked, for a given assignment b, whether or not it is a satisfying assignment), thus showing that relative to a random 



oracle NP is not a subset of subexponential quantum time. In this section, we show that the seemingly small difference 
between the specifications of these two types of queries results in a dramatic change in the query complexity — 0(n 3 ) 
queries suffice to obtain enough information to characterize <£>. Thus black box or oracle techniques do not rule out a 
polynomial time solution to 3SAT by adiabatic quantum search. To reconcile this with the oracle results from it is 
useful to recall that the Cook-Levin theorem, suitably formulated as saying that NP has a iocal-checkability' property, 
does not relativize ^ (see [|l4| for a brief discussion of this issue). In this sense, the results in this section indicate 
that even keeping track about the number of unsatisfied clauses constitutes sufficient structural information about the 
problem to bypass the oracle results. 
More formally, let 

F$(b) := "# unsatisfied clauses in assignment $(6)", 

with b E {0, 1}™. In our black box model, the quantum algorithm is only allowed to access $ via a quantum black-box 
_B$ that reversibly maps \b) |0) i— > \b) \F$(b)). In this section, we prove that the query complexity for 3SAT is 0(n 3 ), 
by showing that F$ is completely determined by its values on the 0(n 3 ) input strings of Hamming weight < 3. Our 
techniques also apply to the Exact Cover problem discussed in [0]. 

For convenience, and without loss of generality, we will not allow repeated variables in the same clause, but instead 
will allow clauses of size less than 3. For example, we can replace the clause {x\ V x\ V X2) with (x\ V X2), and 
(xi V -1X1 V X2) with a constant clause (1) that is always satisfied. Without loss of generality, we can assume that the 
number of such (1) clauses is 0. 

Let us introduce some notation. Let denote the number of clauses in $ that have all three variables without 

negation (e.g. (x\ V Xi V X3)). We will say that these clauses are "of the form" XXX. Let denote the number 

of clauses that have exactly one variable negated (e.g. (x\ V ~^X2 V £3)). Further, we let denote the number of 

clauses that have exactly two variable negated, and denote the number of clauses that have all three variables 

negated. We also define the analogous 1 and 2 variable versions of these expressions. 

Furthermore, if we subscript any of the X with an index, say i, then we only count clauses that have Xi as one 
of the non-negated (or positive) variables. Similarly, if we subscript any of the X with an index, say i, then we only 
count clauses that have ~xl as one of the negated variables. For example, denotes the number of clauses in 

$ that contain the variable Xi and two other positive variables, \XiXjX\ denotes the number of clauses with xi and 
Xj and another positive variable, and denotes the number of clauses that have one of the positive variables 

equal to Xi, another equal to Xj, and another positive variable. The expression \XiXjXk | equals the number of times 
the clause (-1X4 V ~^Xj V Xk) (or equivalent permuted clauses like {xu V ->Xj V ->Xi)) occurs. 

These expressions are symmetric under permutation of the symbols, so for example, \XiXjX\ = \XjXiX\ and 
\XiXjXk\ = \XiXkXj\. 

For example, we have that 

F$(0 tt ) = \XXX\ + \XX\ + \X\ 

since any clause with a negated variable will be satisfied, and the rest will not be satisfied. 
The following definitions will be helpful. For each i E {1, . . . , n} let 

:= \X~XX\ - \XiXX\ + \X~iX\ - \XiX\ + \X~\ - \Xi\. 

For each pair i, j E {1, . . . , n}, i 7^ j, let 

Yij := \X^X\ + \Xa 3 X\-\X~X l X\-\X~X 3 X\ + \X^\ + \X t X J \-\X~X J \-\X~X l \. 

For each triple i, j, k of pairwise distinct integers from {1, . . . , n}, let 

Yijk '■— \XkXiXj\ + \XiXjXk\ + \XjXiXk\ + \XiXjXk\ — \XiXjXk\ — \XjXiXk\ — \XiXkXj\ — \XjXkXi\. 

For each i E {1, . . . , n} let e % denote the string with a 1 in the ith position and 0s elsewhere. For each i,j E 
{1, . . . , n}, i 7^ j, let e y denote the string with a 1 in positions i and j and 0s elsewhere. For each i, j, k £ {1, . . . , n}, 
pairwise distinct, let e % i k denote the string with a 1 in positions i.j and k and 0s elsewhere. 

We now have the next theorem. 



Theorem 2 Let b E {0, 1}™ and let I be the subset of {I, . . . , n} such that b, t = 1 -t=? i E I. Then 

iei i<jei i<j<kei 

Furthermore, 

F»(e') - F*(0") 

F*(e ij ') - F $ (e J ) - F$(e J ) + F*(0 n ) 

F*(e yfc ) - F«(e«) - F $ (e lfe ) - F$ (e jfc ) + F $ (e J ) + F $ (e J ) + F $ (e fe ) - F*(0 n ). 

In other words, in order to be able to evaluate i*$ for every input string {0, 1}™, we only need to query the black-box 
B$ on the 0(n 3 ) inputs with Hamming weight at most 3 (the cases b E {0™, e l , e lJ ', e yfe }). Specifically, we can 
decide whether $ is satisfiable or not by querying the black-box a total of 0(n 3 ) times, after which we use the 
query results to evaluate F$ for all other possible inputs b E {0, 1}™. If any of the strings give F$>(b) = 0, then <f> is 
satisfiable, otherwise it is not satisfiable. (Clearly, with this information we can also answer other decision problems 
like "$ e PP?") The full proof of this theorem is described in the appendix of this article. 

7 Lower Bounds for Adiabatic Algorithms 

In this section we present an easy n-bit problem, for which the adiabatic approach only succeeds if it is allowed 
an exponential delay. We do this by changing an easy problem (the Minimum Hamming Weight Problem) into a 
perturbed version for which the proper solution is as far as possible from its local minimum. It will be shown that for 
this perturbed version, the quantum adiabatic algorithm does indeed require exponential time. 

7.1 The Minimum Hamming Weight Problem 

Consider the adiabatic quantum algorithm that tries to minimize the Hamming weight w(z) of an n bit string z E 
{0, 1}™. We define the initial Hamiltonian by Ho := ^ z w(z) \z)(z\, such that the time-dependent Hamiltonian is 

Hw (t) : = y, w{z)m\ + $ E ™(*M*I- 

ze{o,i}" ze{o,i}" 

As intended, the ground state of the final Hamiltonian is simply |0 • • • 0) with zero energy. 

Since w(z) = z\ H V z n , it is easy to see that H w (t) is a sum of n Hamiltonians, each acting on a single qubit. 

Thus even though H w (t) is a 2" x 2" dimensional matrix, which thus has 2™ eigenstates, these eigenstates and their 
corresponding eigenvalues may be computed by solving the 2 dimensional problem. For the analysis of the minimal 
gap between the two smallest eigenvalues it is again convenient to introduce a relative time-parameter s := ^, which 
ranges from to 1. The eigen-decomposition for the 2 dimensional problem yields: 

- Ao(*)K(*)X«o(*)|+Ai(*)|wi(*)X«i(*)|. (2) 
with 

A (s) = ±-±a/2s 2 -2s + 1 and Ai(s) = ± + ± ^2s 2 - 2s + 1. 

Specifically, at s = we have |u (0)) = |6) = ^-(|0) + |1>) and 1^(0)) = |i) = ^(|0) - |1», while at s = 1 we 

haveMl)) = |0)and|«i(l)) = |l). 

For the n qubit case, it is easily shown that for every y E {0, 1}™ there is an eigenvalue 

X y (s) = (n-w(y))-X {s)+w(y)-X 1 (s), 

where the corresponding eigenvector is the n-fold tensor product 

\v y (s)) := |t> Wl (s))® |%,(s)) ®"-<8> !%,(«)>■ 



Yt = 
Y — 



if i - s s-n 

2 1 s - 1 1 + s J 



Because Ao(s) < Ai(s) for all s, the ground state of H(sT) is |^o( s )> • • • i v o( s )) with eigenvalue tiAq(s). 
The eigenvalues closest to this ground energy are those associated with the w(y) — 1 eigenvectors \v y (s)), which 
have eigenvalue (n — l)Ao(s) + Ai(s). Hence, the energy gap between the two smallest eigenvalues is g(s) = 
\/2s 2 — 2s + 1, with its minimum g m j n = at s = | (i = -j). Because this gap is independent of n, we can 
transform Hq to ii^ adiabatically with a constant delay factor. As a result, the ground state \vq (s)) := |fo(s))® n of 
the system evolves from 1 ■ ■ ■ 0) to 1 ■ ■ ■ 0) in time 0(1). 

We will now discuss an important aspect of the above adiabatic evolution, which we will use in the lower bound of 
the next section. We saw how the initial ground state of the Hamiltonian Hq is the uniform superposition ^= J2 Z \ z ) 
while the final ground state of H w is the zero string |0 n ). Both states share the property that they have an exponentially 
small component in the subspace spanned by computational basis vectors labeled with strings of Hamming weight at 
least (| + e)n. With the eigenvector decomposition of Equation |^ we can see that such an upper bound holds for 
< s < 1. Take for example the vector |1"), which indeed has: 

|<l»|t*00---«o(*)>| < ^r, (3) 

for all s. This bound suggests that a perturbation of the Hamiltonian H w in this subspace will only have an exponen- 
tially small effect on the evolution of the ground state. In the next section we will use this phenomenon to obtain an 
exponential lower bound on the time complexity of a perturbed version of the Minimum Hamming Weight Problem. 



7.2 The Perturbed Hamming Weight Problem 

We will now consider the minimization of a function that is variation of the Hamming weight function of the previous 
section: 

._ / w(z) if w(z) < (| + e)n, 
J[ > '~ \ p(z) if w(z) > (§+e)n, W 

with e > and p(z) a decreasing function that achieves the global minimum f(z) — p(z) = —1 in the w(z) > 
(| + e)n region. Our main result will be the proof that minimum gap of the corresponding adiabatic evolution Hf (t) 
is exponentially small, and hence that the adiabatic minimization of / requires a delay factor that is exponential in the 
input size n. 

For clarity of exposition, we will focus on the special case where 

/w - { w{ i l z z t\;: l {_ (5) 

The proof contains all the ingredients required for the general result mentioned above. 

The fact that this problem is a perturbed version of the Minimum Hamming Weight Problem is best expressed by 

H f (t) := H w (t)-±(n+l)\l")(l«\. 

We will analyze the time-dependent eigenvalues of Hf by comparing them to those of H w . In the previous section, we 
were able to diagonalize the H w matrix by the unitary transformation V(s) that maps the bit string \yi) ® • • • ® \y n ) 
to the tensor product \v yi (s)) (8 • ■ • (8 \v yn (s)). Hence, using s := ^, we have that V^(s) ■ H w (t) ■ V(s) is a diagonal 
matrix with spectrum {X y (s)\y £ {0, 1}™}. By looking at Hf in the eigenbasis of H w we get the following matrix A, 
where we surpress some of the parameters t and s for ease of notation: 

A := ■H w -V-s{n+l)V^\l n ){r\V. 

Note first that for t = and t — T, A is a diagonal matrix. For intermediate values of t, A will have off-diagonal 
entries caused by the perturbation — s(n + 1)|1")(1™| in the Hamiltonian Hf, At t = the minimum eigenvalue is 
zero, which is indicated by the Ax t i = in the top-left corner of the Hamiltonian. At t = T, the minimal eigenvalue 
has changed to —1 (for z = 1"), which coincides with the bottom-right element ^".2" = — 1. The eigenvectors of 
these values are \vq (0)) and|w™(l)), respectively. Intuitively, one expects the critical moment in the time evolution of 
Hf to occur when the ground state has to change from \vq ) to This is indeed the case as we will see next. 



To prove our claim we will introduce another matrix B that equals the matrix A with its entries yl.2,1, . . . , A^n \ 
and A\$, . . . , A\^ erased: 



or, equivalently, 





( M,i 
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B := 





-42,2 


• ■ • A 2 ,2" 
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^2»,2 


• • • A 2 n i2 


* / 


B := A- A\Q n ){Q n 


1 - |o"X 


n \A + 2A 1 


i|0" 



By construction, the state \vq (s)) will be an eigenstate of B for every s with Ai 1 as its eigenvalue. At t = the 
minimum eigenvalue of B coincides with this A\,i = entry; while at the final t = T the minimum eigenvalue (with 
value —1) is 'located' in the (2™ — 1) X (2™ — 1) sub-matrix (corresponding to the subspace orthogonal to |uq (s))). 
Because B transforms continuously between these two extremes, it follows that there is a critical moment s c for which 
the minimum eigenvalue in this subspace and the eigenvalue A\ t \ are identical. In short, at s c the matrix B has a 'zero 
gap' between its two minimum eigenvalues. 

It can also be shown by the definitions of A and V, the fact that V' H W V is diagonal, and the lower bound of 
Equation |]that: 



\\A-B\\ 



< 



V2- \\A\0 n ) 
sV2(n+l) 
s(n + 1) 
V2"" 1 ' 



- (0 n |A|0")(0"|||2 
■\(l n \V(s)\O n )\ 



The optimal matching distance between A and B expresses how close the spectra {Ai, 
of A and B are, and is formally defined by 



,A 2 n} and {fix, . . . ,^ 2 ™} 



d(A, B) 



:= min max |A,- — ivm!, 



with the minimization over all permutations tt G 52« . It is a known result in matrix analysis that for Hermitian 
matrices A and B this distance is upper bounded by \\A — B\\2 (see Section VI. 3 in [||). 

We thus reach the conclusion that for all values of s, the gap g(s) of A (and hence of Hf(s)) will never be bigger 
than the gap of B plus twice the distance \\A — B\\2- At the critical moment s c , when the two minimal eigenvalues 
of B are identical, this implies for the gap of A the upper bound g(s c ) < s c (n + 1)/ V2" -3 , and hence also for the 



Hamiltonian Bf. g n 
delay factor T. 



g 0(-S=). Applying the requirement T 3> g min thus yields the lower bound O(^j-) for the 



7.3 Generalization 

It is not difficult to see that the above lower bound method applies to the larger class of functions mentioned in 
Equation^. The critical property of / is that it only deviates from the Hamming weight function w(z) for those strings 
z that have an exponential small inner-product with the H(s) ground state \vq(s) ■ ■ ■ vq(s)) for all s (the property of 
Equation |]). 

As long as the perturbation p : {0, 1}" — > R in Equation |] is polynomial in n, we have an inequality similar to 
Equation || 

\\(H f -H w )\vo(s)---v (s))h e 2" n ("). (6) 

Hence, if the perturbation p is such that the minimum of / is not /(0"), then the adiabatic algorithm requires a delay 
T g~? that is exponential in the input size of the problem, i.e. T G 2 n ( n \ 



8 Conclusions 



Adiabatic quantum computation is a novel paradigm for the design of quantum algorithms — it is truly quantum in 
the sense that it can be used to speed up searching by a quadratic factor over any classical algorithm. On the question 
of whether this new paradigm may be used to efficiently solve NP-complete problems on a quantum computer — we 
showed that the usual query complexity arguments cannot be used to rule out a polynomial time solution. On the other 
hand, we argue that the adiabatic approach may be thought of as a kind of 'quantum local search' . We designed a family 
of minimization problems that is hard for such local search heuristics, and established an exponential lower bound for 
the adiabatic algorithm for these problems. This provides insights into the limitations of this approach. In an upcoming 
paper [Q], we generalize these techniques to show a similar exponential slowdown for 3SAT. It remains an open 
question whether adiabatic quantum computation can establish an exponential speed-up over traditional computing or 
if there exists a classical algorithm that can simulate the quantum adiabatic process efficiently. 
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A Proof of the Query Complexity Result 

Theorem 2 Let b e {0, l} n and let I be the subset of {I, . . . , n} such that b t = 1 i <E I. Then, 

F(b) = F (0 n ) + ^Y i+ £ Y ij + E 

iei i<jel i<j<k£l 

Furthermore, 

Y t = F(e l )-F(0 n ) 
Y i:j = F(e ij ) - F^) - F{e j ) + F(0 n ) 

Y ijk = F(e ijk ) - F(e ij ) - F(e lk ) - F(e jk ) + F{e l ) + F(e j ) + F(e k ) - F(0 n ). 

Proof: We count the total number of unsatisfied clauses by analyzing each type of clause. 

Firstly, the only clauses of the form XXX that will not be satisfied are those that have all three variables with 
indices in /. This gives us 

\XiXjXk\ 

i<j<k£l 

unsatisfied clauses of the form XXX. Note that if there are less than 3 ones in b then any of the summations over 
i, j, k G I satisfying i < j < k will be empty and thus sum to 0. 

Secondly, the only clauses of the form XXX that will not be satisfied are those that have both of the negated 
variables with indices in I and the positive variable with index not in /. This gives us 

^ \XiXjX\ — ^ (\XiXjXk\ + \XiXkXj\ + \XjXkXi\) 
i<jei i<j<kei 

unsatisfied clauses of the form XXX. 

Thirdly, the only clauses of the form XXX that will not be satisfied will be those that have the negated variable 
with index in I and the positive variables with indices not in /. This gives us 

J2\xxx\- {\x i x j x\ + \x~x i x\)+ (\x~ i x j x k \ + \x~x i x k \ + 1X^X^1) 

iei i<jei i<j<kei 

unsatisfied clauses of the form XXX. 

The only clauses of the form XXX that will not be satisfied are those that contain no variable with index in I. 
This gives us 

\xxx\-J2\x l xx\+ ]T \XiXjX\- ]T \XiXjX k \ 

iei Kjei i<j<kei 

unsatisfied clauses of the form XXX. 
Similarly, we have 

i<jei 

unsatisfied clauses of the form XX, 

iei i<j£l 

unsatisfied clauses of the form XX, 

\xx\-J2\x*x\+ E 1*^1 

iei i<jei 



unsatisfied clauses of the form XX, 



unsatisfied clauses of the form X, and 



E* 

iei 



1*1 -E 1*1 

iei 

unsatisfied clauses of the form X. 

These account for all the unsatisfied clauses. Summing these quantities while rearranging terms according to the 
number of variables in the summations, gives us the first part of the theorem: 

F(b) = \XXX\ + \XX\ + \X\ + Yi + E Y 3 + E Y ^ 

iei i<jei i<j<kei 

= F(o n ) + Y,Y+ E Y v+ E 

iei i<jei i<j<kei 
Notice that for F(e l ) any of the summations with more than one variable will be empty, and we get 

Fie 1 ) = F(O n )+Y t . 
Similarly, for F(e u ) any of the summations with three variables will be empty, and we are left with 

F(e ij ) = FiO^+Yi + Yj+Yij. 

Lastly, for F(e tjk ) we get 

F(e^ k ) - F(0") +Y l + Yj +Y k + Yy + Y lk + Y jk + Y ijk . 
From these equations follow the second part of the theorem. ■ 



